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In this paper we study the steady uniphase and multiphase solutiuons of the 
discretized nonlinear damped wave equation. Conditions for the stability and 
instability of the steady solutions are given; in the instability case the linear 
stable and unstable manifolds are described. 



1. Introduction 

We consider the following mixed problem: 

u tt = {<t{u x )) x + eu xxt , x G (0, 1), t > 0, e ER (1) 

u(0,t) = 0, u(l,t) = P, t > (2) 

u(x, 0) = Uo(x), u t (x, 0) = ui(x), x G (0, 1) (3) 

where P > and a : R — > R is an application with the following properties: 
cr(0) = 0, cr is continously differentiable on R, cr(£) > 0, V£ > 0, cr(£) < 0, 
V£ < 0; there exist a, [3 such that < a < /3 and cr'(^) > 0, V£ e (— oo,a) U 
(j3, oo), cr'(^) < 0, V£ G (a, /3), and a, (3, such that 0<a<f3<f3<j3, such 
that cr(a) = a(/5) and a(a) = <r{0). 

In the paper [1] it is shown that for e — 1 the mixed problem 1-3 has a 
unique solution on the appropiate function space. The uniphase and mul- 
tiphase steady solutions are defined and it is studied the stability of these 
steady solutions. It is presented also a discretization in the spatial variable 
x of the equation (1.1) and some numerical simulations are presented. For 
e = the problem was studied in [2]. 

For cr(£) = a 2 £, a G R we obtain the wave equation studied in a lot of 
papers. 

The aim of this paper is the study of uniphase steady solutions and par- 
tially the study of multiphase steady solutions of the system obtained via 
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a discretization of (1.1) in spatial variable x and then a discretization in 
temporal variable t. 

In the section 2 it is shown that the system obtained from the discretiza- 
tion of (1.1) in the spatial variable x represents an equation Euler-Lagrange 
by rapport with a Lagrange function discrete continuous with damped term. 

The system has a finite number of steady solutions (uniphase and multi- 
phase). 

In the section 3, using the linearized system of the system used in the 
section 2, we study the stability of uniphase steady solutions. We proove 
that such a solution is stable if P G (a, a) or P G {(3,(3) and e > 0. If 
P G (a, (3) the uniphase steady solution is hiperbolic so it is unstable. We 
have determined a polycycle curve for this solution. 

For P G (a, a) or P G ((3,(3), e = 0, we find a curve u c (t), t G R such 
that every of this components is periodic with the period 

7r krr 
T k = — —cosec — 
n^/T In 

For the steady solution 2-phase we present stability conditions depending 
on e, Ti = cr'(a), p = <r'((3), a G (a, a), (3 G {(3,(3). 

In the section 4 it is shown that the discretization of the equation (1.1) 
leads us to a difference system of equations which represents the equations 
Euler-Lagrange with damping for an associated discrete Lagrange function. 
The steady solutions, uniphase and multiphase, are presented and we proove 
that they are in finite number. 

Also, in the section 4 we are studing the stability of the uniphase solu- 
tions, the conditions for a uniphase steady solution to be hiperbolic and the 
linear stable and unstable associated manifolds. 

A similar study for the multiphase steady solutions will be done in a 
future paper. 

2. Semidiscretization of the equation 

U tt = (<t(u x )) x + £U xxt 

For e — 0, the equation (1.1) represents the equation Euler-Lagrange for 
the Lagrange function L : J 2 (R 2 ,R) — > R, L(x,t,u,u x ,u tx ) = \u 2 x — cr(u x ), 
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where J 2 (R 2 , R) is the bundle of the jets of order 2 of the fibrate 7r : R 3 — > 
R 2 . 

Be x k = khi, k = l,2,...,n the division points of the interval [0,1], 
hi — l/n and u k (t) = u(khi,t), k = 1,2, ...,n. For the boundary conditions 
we consider: Uq = u(0,t) = 0, u n = u(l,t) = P. 

The partial derivatives u x , u xx , (cr(u x )) x will be approximate by: 

u x (kh u t) ~ -^-(u k (t) -ttfc-i(t)) 



u xx (kh!, t) ~ ^ (u fc+ i(i) - 2u fc (t) + Ufc_i(t)) (1) 



(<r(wx))a!(A;/ii,t) ~ 

Let be the sequence (u k ,u k ) £ ^R, k = 1, ...,n on the tangent bundle at 
R and L : TR — > R the Lagrange function defined by: 

L(u fc , «*.) = + w (^"( M fc - M fc _i)^) , fc = 1, n - 1 (2) 

where w : R — > R with u>'(0 = The action of L is defined by: 

n 

>l(u,u) = L(uk,u k ) (3) 
fc=i 

In order to obtain the first variation of A we consider the sequence 
(u k (rj), iik(v)) £ TH, k = l,...,n with rj e (—a, a) and Wfc(0) = -Ufc(O) = 
Mfc. The action (2.3) on this sequence is: 



n-1 



k=i 

The first variation of (2.3) is given by 

dA( V ) 



The first variation of (2.3) is: 



= 

r)=0 



u k - ^-a(Au k ) + ^-(r(Aiin) = , k = 1, n - 1 (4) 



where 



and 



Au k+1 = -r-{u k+ i - u k ) , k = l, ...,n - 1 
rii 

Au n = -j-(P - u n _i) 
til 

The system (2.4) represent the semidiscretized system of the equation 

(1.1) for e = 0. 

The semidiscretized system of the equation (1.1) for e ^ is: 

u k -^a(Au k+1 ) + ^a(Au k ) = ^{ii k+1 -2u k + u k ^ 1 ) , fc = l,...,n-l (5) 
til rii tii 

and represents the equation Euler-Lagrange for the function L defined 

(2.2) with the disperssion: 

-To{u k+ i - 2u k + u k ^i) 
tii 

To the system (2.5) we associate the equivalent system on T*R: 

u k = v k , k = l,...,n-l 

1 £ 
ik = ~h~i [ <7 ( AUfe + 1 ) ~~ a (^ U k)} + ^2 ( V k+1 - 2v k + Vk-l) 

with the action: 



m 



(6) 



n-1 



k=l 



+ w(Au n ) 



(7) 



From (2.3) and (2.5) it follows that: 



dA^u " ri 



k=i 



(8) 



From (2.5) and (2.6) it follows that any steady solution of the system 
(2.5) satisfies the conditions: 

a(Au k ) = C, k = l,...,n, J2 TT^k = nP > C e R ( 9 ) 

k=i ^ 



Because the function a is not monotone we obtain two types of steady 
solutions: 

Definition 2.1. A uniphase steady solution u for the system (2.5) is a 
steady solution with the property: Uk = khiP, k — 1, ...,n — 1. 

Definition 2.2. A multiphase steady solution for the system (2.5) is a 
steady solution with the property: u^+i —u^ = hif(k), k — 1, ...,n — 1. 

We denote by E the set of multiphase steady solutions and: 

E + = {ueE : (/(Awfc) > 0, fc = 1, n} 

E~ = {u e E : 3k e {1, n} : <r'(Au fc ) < 0} 
In [1] it is shown that: 

Lemma 2.1. For any steady solution u G E + there is a G (a, a) and 

(3 G (A/9), C G (ct,ct), such that: 

1) ff(Aw fc ) = C, fc = l,...,n 

2) A(« fc ) = a or A(« fc ) = /3, fc = 1, n 

3) ka + (n - k)(5 = nP 

Lemma 2.2. If P < a or P > (3 then E + = 

Lemma 2.3. The system (2.5) has a finite number of steady solutions. 
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3. Properties of the steady solutions of the 

system (2.5) 



The general form of the system (2.6) is: 

F k (u k , u k+ i, iik-i, iik, Uk+i, Uk) =0, k = 1, 2, n 



(1) 



Let u be a steady solution of the system (2.6). The linearized system in 
a neighborhood of u is given by: 



dF k \u + rfwji u + rfwj, U + Yn 



W ; 



j'=i 



where 77^- G (—a, a), j = l,...,n. 

The linearized system associated to (2.5) is: 



= 



j = l,...,n 



hi 



[w j+1 - 2wj + Wj-i), j = l,...,n-l 



A solution for the system (3.3) is: 



7r/c 



n 



(2) 



™3 - ^(Au j+ l)(%l - Wj) + ^<T'(AUj)K- - Wj-l) = (3) 



lUj = exp (Xj)exp (ia k j), a k = — , k = 1, ...,n - 1, A e C (4) 



From (3.3) and (3.5) yields: 

\ 2 exp ia k - -^a'(Au k+1 )exp ia k (exp ia k -i) + -La'(Au k )(exp ia k -i) 



£ 

hi 



2 X(exp ia k ^i) 2 



(5) 
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Theorem 3.1. If P G (a, a) or P e /3) the uniphase steady solution 
Uk = khiP, k = l,...,n — 1 of the system (2.5) is asymptotically stable if 
e > and unstable if e < 0. 

Proof: For the uniphase steady solution -u, the equations (3.5) become: 



1 £ 

A 2 exp iak — T2 T ( ex P i a k-i) 2 — joK^P i a k-i) 2 = 0, k — 1, n — 1 (6) 



where r = cr'(P). We denote: 

— : 

exp tdk 

The equations (3.6) become: 



a ■ 2 T^k 

= —Asm — 
2n 



A 2 — e/ifcn 2 A — /ifen 2 r = 0, fc = 1, n — 1 



with the roots: 



1 



\ 2 x 2 

A = 



e ± e + 



4t 
p fc n 2 



, fc = 1, n — 1 



(7) 



(8) 



If £ > 0, because p^ < yields that Re A& < so it follows that the 
steady uniphase solution u is asymptotically stable. 

If e < 0, it follows that Re \ k < so the steady uniphase solution is 
unstable. 

Theorem 3.2. If P e (a, /?) the uniphase steady solution -u of the system 
(2.5) is hyperbolic so it is unstable. 

Proof: For the uniphase steady solution u in the equations (3.7) become: 



A 2 — e/ifcn 2 A + p&n 2 p = 0, k — 1, n — 1 
where p = — cr'(P); the solutions of these equations are: 



A fc = \vkn 2 



(9) 



(10) 



a; 



:Pfc™ 



V Vkn A 



, /c = 1, n — 1 
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The real solutions A + , A~ have the property: A + ■ A~ = p k n 2 p < 0. It 
follows that they have opposite signs. The steady solution is unstable, it is 
hyperbolic. 

The linear manifolds associated to the steady solution of hyperbolic type 

are: 

n-l 

u k (t) = h\kP + exp (t\ k ) ^ r/jexp (ia k j), t e R, k — 1, n — 1 (11) 

i=i 

n-l 

■u^(t) = /ii/cP + exp (t\ k ) r\jtxp (ia k j), t e R, = 1, n — 1 (12) 

j'=i 

and have the properties: 

lim Mfc(t) = u k , lim ut(i) = Uk, k — 1, n — 1 

t— >— oo i^oo 

Theorem 3.3. The curves ujyfc +1 , A; = 1, 2, n, t G R given by: 

*%(*) u u k+1 (t ) 

are policycles for the uniphase steady solution, that is 

hm u h kk+1 (t) = u k 

Proof: From (3.12) and (3.13) we have: 

u k u u (t) 

lim u h kk+1 (t) = - u k+l + u k lim 

From (3.11) yields: 

n 

hi(k+ l)P + exp(t\t +1 ) J2 Vjexp(ia k+1 j) 

lim = lim ^ 

fc hi(k)P + exp(t\ k ) ^2 Vj ex P(^ a kj) 

3=1 

n 

Y,Vjexp{ia k+1 j) 

Ai.i 



Mfcfc+l - — — 77\ Uk +i + r e K (Id) 



^hme*p(*(A+ +1 -A+))^ 



J2vjexp(ia k j) 

3=1 
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From (3.10) we have A^ +1 < X k and A fc+1 < \ k . It follows that 

■ u h k {t) 



lim%^ = 



Thus we obtain: 

lim u h k+1 {t) = u k 

t — >oo 

Using (3.11) and (3.13) it follows: 



hrn^u^t) = u k+1 + u k+l t hm ^-^ 
From (3.12) yields: 



J2vjexp(ia k j) 
j=i 

n 

J2vjexp(ia k+1 j] 

3=1 



Thus we obtain: 



t ]im o M^ +1 = u k+l 
The curves (3.13) are a polycycle considering: 

■U n+ 1 = Mi 

Theorem 3.4. If P e (a, a) or P e /3) and e = the curves: 



Ufe(i) = Mfc + e:rp (A c t) rjjexp (ia k j), k — 1, ...,n, t G R (14) 

7T A/7T 

are periodic, with the period = — -=sec — , k = l,...,n, where A c = 

n^/T 2n 

kn 

±2inJrsin — . 

2n 



Proof: For e = the equation (3.7) is: 

A 2 — fikn 2 T = 0, k — 1, n — 1 

with the roots: 

A c = ±2in\frsin — 
2n 

The uniphase steady solution uu = kh\P is center for the system (2.5) 
with e = 0. From the condition u%(t) = u% +1 (t + 7^), k — 1, ...,n — 1, it 

A/7T 

follows that errp (A c (t + 7^)) = errp (A c (t). Hence: 2riy / rsin — Tk = 2n or 

Zilh 

7T A;vr 
T fc = — — cosec— 
nyr In 

We continue our study with a multiphase steady solution u G To be 
more specific we will consider only 2-phase steady solutions. From Lemma 
2.1. it follows that for a G (a, a), (3 G ((3,(3), with 7(a) = {i G {1,2, ...,n}, 
A-Uj = a}, I((3) = {i G {1, n}, A-Uj = /3} we have 

7(a)U7(/3) = {l,...,n} 

carrf 7(a) • a + (n — carci I{(3))(3 = nP 
Let consider the following sets: 

7(a) = {iG 7(a) : i + le 7(a)} 

J{(3) = {i G I{(3) : 

and 

7T(a,/3) = {i G 7(a), i + le 7(/3)} 
From (3.5) we obtain: 

A 2 + Aen 2 sin 2 ^\ + 4n 2 sm 2 ^r = 0, p G 7(a) (15) 

A 2 + Aen 2 \sin 2< ^ + 4pn 2 W^r = 0, p G J(/3) (16) 
\ 2 +Aen 2 \sin 2 -^-—n 2 [p(exp ia p ^i)+r{exp {—ia p ) — l\ = 0, p G K(a,/3) (17) 
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where r = cr'(a), p = <t'((3), p > r. 
Lemma 3.1. 

a) . The equations (3.15) have complex roots with negative real part if 
and only if: 

e > and e < q-, r = max J (a) (18) 

nsin— 
2 

The equations (3.15) have real negative roots if and only if: 

e > and e > q— , q = min J{a) (19) 

nsin— 
2 

b) . The equations (3.16) have complex roots with negative real part if 
and only if: 

e > and e < v ^ Qr , r = max J(/3) (20) 

nsin— 
2 

The equations (3.16) have real negative roots if and only if: 

e > and e > ^ a , q = min J{f3) (21) 

nsin— 
2 

c) . The equations (3.17) have complex roots with negative real part if 
and only if: 

P + T-^fK \jp + T+^K 



z > Vr , a q < ^ < V , a r (22) 

nsin— nsin — 

2 2 

q = min K(a,j3)), r = max K(a,/3)) 



where 

A = (p + r) 2 -(p-r)W| 
Proof: From the equations (3.15) and (3.16) it follows directly a), and 

b). 

The equations (3.17) may be written in the following form: 

\ 2 +Aen 2 sin 2 ^\+2n 2 (p+T)sin 2< ^--in 2 (p-T)sin a p = 0,p e K(a,/3) (23) 
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The necessary and sufficient condition for these equations to admit com- 
plex roots with negative real part is: 



8 I6e 2 n*sin^ 

2 



2n 2 (p + r)sin 2 ^- > 



(24) 



and 



Aen 2 sin 2 ^- > 0, 2n 2 (p + r)m 2 ^ > 



Theorem 3.5. The 2-phase steady solution is stable if one of the follow- 
ing conditions holds: 

a)- . 

Jp + r-^/A y/T . 
q < e < , q = mm K (a, p)) 



b). 



nsm— 
2 



n 



< e < q , r = max K (a, p)) 



nsm- 



7T 

'2n 



nsm- 



max 



P 



P + t 



nsm 



a q2 ' 



nsm 



"qii 



A 1 . 

< e < mm 



P + T 



On ' 



nsm- 



2 2 / \ 2 

g 2 = rain J(/3), q 3 = min K(a,j3)) 
T\ = max J (a), r 3 = max K(a,j3)) 



nsm- 



d). 



max 



p + r-VK 



nsm 



a qi > 



nsm 



"93 



< £ < mm 



p + r + 



2 2 / v 2 

<7i = min J (a), q^ = min K(a,(3)) 



nsm- 



nsm- 



a r , 



J 
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r 2 = max J(/3), r 3 = max K(a,j3)) 

The proof is a direct consequence of the Lemma 3.1. 
Theorem 3.6. If e = the curves 

n 

M pW =u p + exp (A^t) Vkexp (ia p k), p G J (a) 

k=l 
n 

Up(t) =u p + exp (A^) Y Vkexp {ia p k), p G J{(3) 

k=l 

where A£, p G J(a) respectively A£, p G J(/3) are roots of the equations: 

A 2 + 4n 2 sin 2< ^-T = 0, p G J (a) 

A 2 + 4n 2 sin 2 % = 0, p G J(/3) 
are periodic with the period 

T p = — -=cosec—, p G J (a) 



respectively 

7T p7l 

T p = — -cosec—, p G J(f3) 
The proof follows from the above theorem. 

Remark. For e < the equations (3.15) and (3.16) have positive roots 
and it follows that the 2-phase steady solution is unstable. 

For e < or e not satisfying one of the above conditions a)., b)., c)., the 
roots of the equations (4.10), (4.11), (4.12) may have positive or negative 
real part in the case they are complex, or they are positive or negative real 
numbers. In this case the steady 2-phase solution is unstable; it is hyperbolic 
saddle. The associated linear manifolds may be described completely in the 
same manner we have done for the uniphase steady solutions. 

4. The discretization of the equation 

u u = (cr(u x )) x + eu xxt 
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Let x k = khi, k = l,...,n the division points of the interval [0,1] with 
hi = 1/n and t p = ph 2 , p — 1, ...,m the division points of the interval [0, 1], 
h 2 = 1/n. 

We'll approximate the derivates Uk(t), v>k(t) by — [u k — ufT 1 ] and T^\u p k +1 — 

h 2 h 2 

2u k + Ufc -1 ], where u p = u(kh\,ph 2 ). 

We call the discrete system associated to the equation (1.1) the following 
system: 

1 r p+l ~ »> n-li 1 /I 

■u 



ft? 1 "* 



fc = 1, n — 1, p = 1, to — 1 
The corresponding boundary conditions are: 

4 = 0, < = P>0, p = l,...,m (2) 

The system (4.1) represent the discrete Euler-Lagrange equations for the 
discrete Lagrange function: 

L(uU,ul,ul +1 ) = ^K +1 -ulf + w(±(ul-uU)) (3) 

where w : R — > R, = cr(£), and the dispersive term is: 

[uSi ~ ~ 2 ( M fc +1 - «fc) + - U k-l] 



h\h 2 

For (4.1) it follows that the steady solutions u = u k satisfy: 

1 _ 1 

^(^■K+i-Wfc)) = <7 (^-K- u fc-i))> fc = l,-,n-l (4) 



From (4.2) we obtain: 



14 



Thus a steady solution u satisfies the following: 

u\ - uU = P- uU, a(±-(u p k - u p _J) = C, k = 1, ...,n, CeR (5) 

Definition 4.1. A steady solution u = (u k ) with u p k = khiP, k = 
1, n — 1, p = 1, m is called uniphase steady solution. 

A steady solution u = (u k ) with the property u p k — u k _i = hif(k), k = 
1, ...,n — 1 is called multiphase steady solution. 

Let E = {u = (u p k )} the set of multiphase steady solution. We consider 
the following sets: 

E + = [u = (u p ) : a'(±-(u p k - uU)) > 0, k = 1, n - l} 

E- = [u = (u p ) : a\±Au p k - u^)) > 0, k = 1, n - l} 

Theorem 4.1. The system (4.1) has a finite number of multiphase steady 
solutions. 

n 

Let u = (u k ) a steady solution of (4.1) and u k = u k + ^^wf(A;), /c = 

2=1 

1, n — 1 the components of a vector with 7/; G (—a, a). 

The linearized system associated to the system (4.1) in a neighborghood 
of the steady solution u is: 

hl(wf +1 - 2wf + wf' 1 ) - h\h\ p(^K+i - ul)(wf +1 - wf)- 

^(^("fc - u k-i)( w f ~ w f-i) 

= eh 2 [w p +l + wf +1 - 2« +1 - wf) + wtl - wU (6) 

p — 1, m, I — 1, n — 1 
A solution of (4.6) has the form: 

wf = X p exp (idkl), a,k = — , k = l,...,n, A G C (7) 
Replaing in (4.6) we obtain: 
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where 



h\(\ 2 — 2A + l)exp (ia,k) — h\h\(exp 2iak — exp ia k ) — 
— s k \(exp ia k — 1) — eh 2 (\ 2 — A) (exp ia^ — l) 2 = (8) 
k — 1, n 

From (4.8) it follows: 

y?[h\ — eh 2 (exp ia k — l) 2 ] — \[2h\exp ia k + h\hi(s k +\exp ia k )(exp ia k — 1) — 

— Sk(exp iak — 1) — h 2 e(exp iak — l) 2 ] + h\exp iak = (9) 
k — 1, n 

For a uniphase steady solution u = (u k ), u k = kh\P and r = cr'(P) > 0, 
the equations (4.9) become: 

\ 2 (h\ - eh 2 fi k ) - \(2hj - f%h X Tn h - eh 2 /i k ) + h\ = (10) 

k — 1, n 

where 

(exp (ia k - l)) 2 A ■ 2 k7C 
exp — 1) In 

h 2 

We suppose for the rest of this section that hi,h 2 € (0, 1) and e ^ 1 



h 2 fik 

Theorem 4.1. The uniphase steady solution for the system (4.1) is as- 
simptotically stable if and only if the following condition hold: 

Proof: The equations (4.10) have the roots in modulus less than 1 if and 
only if (4.11) hold. 

A necessary and sufficient condition for the k-th equation in (4.10) to 
have the modulus of its roots less than one is: 

h 2 

hl 1<0 



h\ - eh 2 [i k 
16 



_ , _ hi - eh 2 Hk \2hl + h\hiT^ k - £h 2 /j, k \ h\ - eh 2 ^ k 

h 2 2 < h\ < + h\ 

The inequatities (4.12) hold if and only if 



4/ii / 2/i? /ii/i 2 \ 

But /ife > so it follows that (4.13) hold for any A; = l,...,n if and 

only if (4.11) hold. 

Theorem 4.2. The uniphase steady solution u = (u p k ) of the system (4.1) 
is unstable if one of the following conditions hold: 

a). 



hi f h{ h\h 



2 



b). 



C). 



T< h\sin^ £e \ °°' 4W^ + 2 T , 



fei / fej /if N 



Proof : a). The equations (4.10) have the roots in modulus greater than 
one only if a), hold. A necessary and sufficient condition for the k-th equation 
from (4.9) to have roots in modulus greater than one is: 

hj - eh 2 n k _ 
h\ 



hl-eh 2 fi k \2h\ + hl^rfik - eh 2 fi k \ h\ - eh 2 fi k 

hi hi hi + 1 ' 

From (4.14) we obtain 

4/ii / 2h\ h x h 2 \ 

h 2 n k V h 2 jj k 2 
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b)., c). can be prooved in an analogous manner. 

Theorem 4.3. The uniphase steady solution u = (u k ) for the system 
(4.1) with e = is unstable if 

T < 



h\sin 2 ^ 

Proof : For 5 = 0, the equation (4.10) become: 

h 2 

\ 2 -\(2 + -^rfi k ) + l =0, k = l,...,n (16) 
hi 

h 2 

From (4.16) we have AiA 2 = 1 and Ai + A 2 = 2 + -r^T[i k . A necessary and 

hi 

sufficient condition for k-ih equation from (4.16) to have one root greater 
than one and the other less than one is: 

T < 



hnsin 2 ^- 



Proposition 4.1. Let u = u p k a uniphase steady solution for the system 
(4.1) which satisfies one of the condition b). or c). from the theorem 4.2. 
and A^, \ k the roots of the k-th. equation (4.9) with \\ k \ > 1 and |A^| < 1. 

The linear stable manifold is given by: 



u 



k = kh iP + (A fc f Vi^xp (ia k l) 



p—l,...,m, k — l,...,n — l 
The linear unstable manifold is given by: 

n 

u 7 = kh i p + (KTj^Viexp (ia k l) 
i=i 

p—l,...,m, k = l,...,n — l 
In these conditions it follows that: 

lim ut p = khiP and lim ut p = khiP 
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